/*=========================== = =====*/ 

/*= = = = = = = = = = Compute Active Periods = = = = = = = = = = =*/ 

/* = = = = = = = = for an Elementary Video Stream = = = = = = = = =*/ 

/* = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = */ 

/* Index / denotes the i-th frame, q\ denotes the queue size at the instant 
just after the i-th frame arrival and index j denotes the ;'-th active period. */ 

1 /* Perform initialization */ 
1.1 $ 1 r <-l//,g 1 «-</ 1 ;i<-y<-l; 

2 

2.1 If ( [(q,)l( r)] < Vf) /* no backlog present at the start of next frame */ 

2.1.1 q i+ \ <r- d i+ \; 

2. 1 .2 tj <— Vf + [(qi)l( r)]; /* compute end time of current active period */ 
2.1.3 

2.1.4 sj r <- I* start time of new active period */ 

2.2 else 

2.2. 1 q i+ \ <r- q t + d i+ \ - Vf, /* there is backlog carried to current frame */. 

2.3 endif 
3 

3.1 

3.2 If 0 < N) 
3.2.1 goto Step 2; 

3.3 endif 
4 

4.1 s(r) <— maxi < £ < ^q^; I* compute maximum queue length observed */ 

4.2 tj <— N /f + [(<7#)/( r)\, I* compute end time of last active period *l 

4.3 n a (r)*- j; I* store the number of active periods */ 

4.4 Define s„ fl(r ) + i = oo; 
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/* INIT. B holds the available buffer, L B holds the sum of the unrestricted queue sizes of 
the busy 

periods with loss at the time instants at which the maximum queue occurs, «/ holds the number 
of busy periods with loss, % the number of busy periods, and n a the number of active periods. 

1.1 5<~oo 0;n/<— 0;M<— stream_size; 

/* Compute busy periods for rate P, and buffer size B. */ 

1.2 process_activejperiods(l, n a , B)\ 

/* Insert into the heap the buffer point at which a loss occurs 

for each busy period, compute the buffer size for the first break for 

each busy period and insert it into the heap as well. */ 
13 For b = 1 to rib /* f° r eac h busy period b */ 

1.3.1 heap_insert(0£*, LOSS, b)\ 

1.3.2 5'<-compute_next_break(Z>); 

1 .3.3 heap _insert(5', BREAK, b)\ 
1.4 endfor 

2 /* Extract the maximum buffer from the heap, and process 
the corresponding busy period until the heap becomes empty. */ 

2.1 (b, B, cause) <- heapextractjnaxO; 

/* Let (fy, t p ), (Sp+\, tp+\), (s qi t q ) be the active periods contained within the busy period b. 
Also, let (s py t p ), (sjj tj) be the active periods contained within the interval ( a^, T£>). */ 

2.2 If (causes LOSS) 

/* update n\ and Lg variables */ 

2.2.1 L B <^L B + ( If = y r(ti *j) - P(tj - s p ) ); 

2.2.2 «/<— w/+l; 

2.3 else /* cause = BREAK */ 
/* update n\ and L$ variables */ 

2.3.1 L B ^L B -{ 1/ = y r{t t -si) - P(tj - s p ) ); 

2.3.2 «/^w/-l; 

/* Process the break in busy period b: Compute the new busy periods, and update n\ and L B 
when a busy period already experiences loss. For the remaining new busy periods, insert 
into the heap the buffer points at which the first losses occur. Finally, for all the new 
busy periods, compute the buffer sizes for the first break and insert them into the heap. */ 

2.3.3 process J>reak(6, 5); 

2.4 endif 

2.5 output_point(5,[(L 5 -n/5)/(A^)]); 

2.6 If (heap not empty) 
2.6.1 goto Step 2; 

2.7 endif 

/* Output the last point of the loss curve */ 

2.8 Output (0, [(r-P)/(r)]); 

2.9 STOP; 
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